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ABSTRACT 


Akima cubic spline and spatial Fourier series (SFS) techniques for modeling ocean 
bottom contours from bathymetric data were comparatively analyzed. SFS methods en- 
countered difficulty in ocean bottom reconstruction despite several enhancements to the 
fundamental technique. Both Akima cubic spline and SFS approaches proved unsatis- 
factory in reproducing first and second-order derivatives for several arbitrarily shaped 
ocean bottom contours. 

A simple reflection angle algorithm for arbitrary one-dimensional bottom models 
was examined and found to be accurate. A graphical demonstration of acoustic ray 
interaction with the ocean bottom using a variety of mathematical ocean bottom models 


and the reflection angle algorithm: was conducted. 
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I INTRODUCTION 


The primary endeavour of this thesis research was the comparative evaluation of 
existing numerical cubic spline techniques with spatial Fourier series (SFS) methods for 
the interpolative modeling of real, discrete, ocean bathymetric data. The Akima cubic 
spline was chosen for its reported ability to produce ocean bottom contour curves which 
pass through the original data while minimizing artificial oscillations in the interpolation. 
The Akima cubic spline logic was implemented via the IMSL Version 10 double preci- 
sion computer program DCSAKM. First and second-order derivatives of the resulting 
spline were calculated from a second IMSL Version 10 double precision program 
DCSDER. The SFS technique was developed as a separate generic computer program 
subroutine. 

Seven differentiable inathematical bathymetric data models were employed in deter- 
mining the accuracy and relative worth of both the Akima cubic spline and SFS repres- 
entations of the ocean bottom contour, including the first and second-order derivatives 
of these contours. A series of numerical enhancements were also applied to both mod- 
eling techniques to optimize their utility in accurately reconstructing sampled ocean 
bottoins. 

The second objective of this research was the development and evaluation of a re- 
flection angle algorithm, derived from the principles of the calculus as outlined by 
Goodman [Ref. 1], which enables the computation of reflection angles for acoustic rays 
incident upon an ocean bottom whose variation in depth is modeled as a function of 
down range. This algorithm requires exact first derivatives for angle calculations. Hence 
the apriori requirement for accurate ocean bottom modeling. 

Both of the developments outlined above were then incorporated into an existing 


acoustic ray trace program aud several simple case studies investigated. 


Ii. THEORETICAL ANALYSIS 


A. SPATIAL FOURIER SERIES REPRESENTATION OF OCEAN BOTTOM 
CONTOURS 

For any ocean region of horizontal range extent Z meters with a one-dimensional, 
range dependent bathymetry, the ocean bottom contour y,(z) can be modeled by the 


following spatial Fourier series (SES): 


S 
ув(2) = 1001 + 2» Iw,| cos(2nsfzz — £w)) ,O&z & Z. 29 


mz [| 


The magnitude and phase angle of the complex Fourier series coefficient w,, for spatial 
harmonic s, are given by |w,| and Zw, , respectively. The average value for depth, which 
corresponds to zero spatial frequency or the “DC component” of the bathymetry data, 
is represented by |]. The fundamental spatial frequency f, is defined as: 


1 
fa "NT; 





(2.2) 


where NN is the total number of bathymetric data points and T, is the uniform sampling 
period in the z direction in meters per sample. The highest spatial harmonic is S and 1s 
governed by the number of data points: 


es = , for even N (2°) 


and 


N— 1 


=== 





, for odd N. (2.4) 


The first and second-order derivatives of the SFS representation of discrete ocean 


bathymetry are readily determined from Equation (2.1): 


S 
ya C) = —4nfy,> s|w,| sin(2zsfxz — £w) ,0< z< Z (2.5) 


=|| 


and 


S 


Ув (2) = — а, | и, | : ол θε Ζ. (216) 


s5 


The complex Fourier series coefficients w, are determined from the spatial discrete 


Fourier transform (DFT) of the ocean bathymetric data as follows: 


N—1 
να Eu ees 
wm DPTR, s ο (2.7) 


The term H7 is known as the DFT weighting factor and is given by 
Еа а (2.8) 
here 
№= 25 + 1 τυ) 


when the ocean bottom sampling frequency equals the Nyquist rate. 

The reconstruction of anv one-dimensional function or contour by a Fourier series 
of finite extent will be characterized by small oscillatory deviations from the exact 
function values in the interpolated region between the discrete samples. These oscil- 
lations are often of little consequence due to the small amplitudes involved. However, 
efforts to model a simple ocean bottom profile of constant non zero slope by SI'S 
methods generate significant amplitude errors at the function extremeties. This develop- 
ment highlights a particular problem that SES methods have with apparent function 
discontinuities represented by the boundaries of the discrete data string. This anomoly, 
first documented by Willard Gibbs in 1899 and known as "Gibbs. phenomenon" [Ref. 
2], can be partially redressed by the inclusion of more data samples. This results in the 
compression of the larger amplitude oscillations toward the discontinuities. However, the 
peak amplitude of these ripples, being independent of the number of discrete data, re- 
mains constant. The problems associated with SIS methods and the "Gibbs! phenome- 
non” take on greater significance when derivatives of the bottom contour are sought. 
The ripples generated at the extremities by the Fourier series representation will be 


magnified in subsequent calculauons making reliable deternunation of function 


derivatives difficult. Strum and. Kirk. (Ref. 2] provide a more detailed treatment of 
“Gibbs’ phenomenon’. 

Lanczos [Ref. 3] argued that the impact of this phenomenon could be minimized by 
hastening the convergence of the Fourier series and proposed the following modified 


complex Fourier series coefficient w,, as a means of achieving this: 


sin(sz/ S) 
ваа 


п = IS м; , 5 5 0 (2.10) 


where w, 1s the original coefficient given by Equation (2.7). 

An additional modification to Lanczos’ formula was suggested by Hamming [Ref. 
4| and mvolved increasing the harmonic S by lI. The final expression for spatial complex 
Fourier series coefficients emploving the Lanczos smoothing methodologv therefore be- 
comes 

n(sz/(S 4- 1)) 
Ws = aS) > we, SFO. emer 
B. THE REFLECTION ANGLE ALGORITHM 

Acoustic rays obey the reflection laws applied in optics. Equation (2.12) relates the 

angle of incidence € and angle of reflection 2’ for a sound ray innpinging the ocean sur- 


face or ocean bottom: 
C=C. σα 


Consider a rav launched from a sound source with initial angle f, and propagate 
in an isospeed ocean with a perfectly flat sea surface. When the ocean bottom is flat, the 
angle of incidence { will remain constant for successive reflections as shown in Figure 


l. The magnitude of the reflection angle 15 given by 
S — fig. (2.13) 


When the ocean bottom shoals (up slope) with constant angle of slope y, the re- 
lationship given by Equation (2.13) is invalidated and subsequent acoustic ray reflections 
become linearly dependent upon the bottom slope. This observation, reported by Officer 


[Ref. 5] and illustrated in Figure 2, is given by 


Cn = ο ο) ο ο (2.14) 


where (η = 1) corresponds to the first bottom bounce with higher values for m repres- 
enting subsequent surface and bottom reflections. Negative values for ¢, indicate acous- 
tic rays which have been reflected back toward their source. Asstgning a negatuve value 
to y in Equation (2.14) yields the relationship for §, after an initial reflection from a 
bottom contour that deepens (down slope) with constant angle of slope y. Figure 3 aud 


Equation (2.15) refer to this situation: 


om tinct DN n= 1.23... (2.15) 
κ Surface 
| 
| 
s | ὃν 
| 
| Bo | 
| | 
| 
A Bottom 
Where C=C =By 


Figure ]. Acoustic Ray Impinging a Flat Ocean Bottom: This diagram represents 
an acoustic ray launched with mitral angle of propagation f, in an 
isospeed ocean with a flat ocean surface. 


Surface 


Bottom 





Figure 2. Acoustic Ray Inpinging a Constant Up-Slope Ocean Bottom: This dia- 
gram represents an acoustic ray launched with initial angle of propa- 
gallon Û, in an isospeed ocean with a flat ocean surface. 
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Figure 3. Acoustic Ray Impinging a Constant Down-Slope Ocean Bottom: This 
diagram represents an acoustic rav launched with inital angle of propa- 
gation fin an isospeed ocean with a flat ocean surface. 


These laws of reflection are equally valid in more general ocean models. However, 
their application is correspondingly more complex. In the case of constant sound speed 
c, acoustic rays travel in straight lines within the initial YZ launch plane making the 
calculation of £ relatively simple. However, if sound speed is made to vary with depth 
such that 


С = c(y), (2.16) 


then the acoustic ray will remain in the initial YZ launch plane but travel a curved path. 
Additionally, if the ocean bottom contour y, varies arbitrarily with down range z from 


the sound source where 
YB = ув(2), (2.17) 


then the relationships given by Equations (2.14) and (2.15) are no longer relevant since 
the angle of slope y does not remiain constant. 

The calculation of € in an ocean model characterized by Equations (2.16) and (2.17) 
is achieved by representing the direction of the acoustic ray at any pomit along its path 


by the unit vector 


η... με (2.19) 


as shown in Figure 4. 


Acoustic Ray 





Figure 4. Vector Representation of an Acoustic Ray: Ihis diagram shows the 
vector representation of a sound ray travelling in an ocean medium 
where sound speed is a function of depth y. 


The vector components :(y) and w(j) are the direction cosines with respect to the Y and 


Z axes, respectively, describing the direction of propagation of the ray and are given by 
v(y) 2 cos f(y) (2319) 
and 
w(y) = sin B(y). (2.20) 


In addition, the tangent plane and normal vector at any point P(x,,z,) on the ocean 
bottom must be determined. 
The tangent plane at a particular point Р(х,у.) ОП а two-dimensional surface 


given by 
ΡΞ (2.21) 


can be determined as follows (see Goodman [Ref. 1]): 


στοκ aa) zO ee (222)‏ ر 
where‏ 
fle) = EI (2.23)‏ 
χνϐ»-ὁ ôx τ‏ 
and‏ 
Vno m ten (2.24)‏ 
ez (1 1 xo Xo Z2 Zg'‏ 210.40 


Goodman further showed the unit vector normal to the tangent plane al xX o aE 
be 


— 





^ My 
П = (2.25) 
ln, | 
where 
ny =—f,(xg,29)i +7 —L(xg.29)k (2.26) 


and 
OECD IO. (2.27) 


Substitution of the function y, given by Equation (2.17) into Equation (2.22) yields the 
following general expression for the tangent plane at any point P(j%,z,) on an ocean 


bottom which is an arbitrary function of down range z. 
J — yo = Y'a(z)(2 — zo). (2.28) 


Similarly, substitution of the function y, into Equation (2.25) provides the following 


general relationship for the vector normal to the bottom at any point P(y,,z,): 


— 


R n 
ر‎ = (2.29) 
where 
Hg, — +j — y B(Zo)k (2.30) 
and 
ngil — A 1 0" в(20)) . (2.31) 


The angle of slope y at any point along the one-dimensional ocean floor y, can be 


determined from elementary vector calculus with the use of Equation (2.29) as follows: 
ATE 
y=cos (yep) (2:37) 


where y is the unit vector lying along the y axis. The vector representation of an acoustic 
ray given by Equation (2.18) can be used in defining the sound ray n, incident upon the 


ocean bottom where 
A A A 
n; = v, J + м; К. (2533) 


The angle of arrival f, which is measured from the y axis to n, and is necessary in the 


calculation of the angle of incidence £, is given by 


f; —cos (p.n). (2.34) 


Finally, with. general expressions determined for y and fj, the angle of incidence C 


and the launch angle for the reflected ray fj, can be evaluated from the following re- 


lationships: 
6 — fi; — Aygy (2.35) 
and 
В. = 180° – (6 + Дуру) (2.36) 
where 
ور‎ = 220. (2.37) 
ly” (za) | 


Figures 5 and 6 illustrate the respective cases of an acoustic ray incident upon an ocean 


bottom with negative and positive gradients, respectively. 
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Figure 5. Acoustic Ray Impinging an Arbitrary Ocean Bottom: This diagram re- 
presents an acoustic ray incident upon an ocean bottom which 15 
shoaling (up slope) at the point of impact. 
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Figure 6. Acoustic Ray Impinging an Arbitrary Ocean Bottom: This diagram re- 
presents an acoustic ray incident upon an ocean bottom which 15 
deepening (down slope) at the point of impact. 


lil. OCEAN BOTTOM MODELING OF DISCRETE BATHYMETRY 


A. A COMPARATIVE ANALYSIS OF AKIMA CUBIC SPLINE AND SPATIAL 
FOURIER SERIES REPRESENTATIONS OF BATHYMETRIC DATA 

The Akima Cubic Spline (ACS) algorithm accurately models most continuous 
functions from discrete data but will frequently return inaccurate first and second-order 
derivatives. Reliable derivative calculations are crucial for reflection angle and ray 
acoustic computations and consequently an alternative modeling technique was needed. 
One of the primary focuses of this thesis, therefore, was the investigation of Spatial 
Fourier Series (SFS) methods for modeling the ocean bottom from discrete bathymetric 
data. This particular approach possessed the added attraction in that bottom roughness 
could be related to spatial frequency components. Furthermore, derivative calculations 
could be readily determined from the resultant series representation of the contour ob- 
viating the requirement to call complex external computer programs. 

Synthetic bathymetric profiles comprising 1] equidistant data points over a 1000 
meter range were generated from seven exact mathematical functions containing features 
consistent with real ocean bottoms. ACS and SFS representations of each bottom con- 
tour, including their first and second-order derivatives, Were then computed. Percentage 
error calculations were made with respect to the exact function values for each model 
and its associated derivatives. When the difference between the modeled and exact data 
was of magnitude less than 10°°, the percentage error was assuined to be zero. Similarly, 
percentage errors greater than 106 were represented in tabular output as a series of as- 
terisks. 

The exact functions used to generate the artificial ocean bottom contours were: 

e Flat bottom ( v,(z) 2 100 m ) 

* Constant downslope bottom ( },(z) = 220.0+ 0.22 m_ ) 

e Constant upslope bottom ( y,(z) = 420.0 —0.2z m ) 

e One period cosine bottom ( },(z) = 450.0 — 250.0 cos(0.002zz) т ) 


i ; . 125000.0 
УЙ — ). Á" ڪا --- - و کس‎ OG Mě 
e Modified Witch of Agnesi bottom ( yg(z) = 100.0 — (0.012? + 2500.0) m) 


e Half catenary bottom ( y,(z) = 315.0 — 30.0(e7/68-9 + @-2/600.0) m ) 


(z — 5000.0) 


¢ Parabolic bottom ( y,;(z) = 50.0 + 2500.0 
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Figures 7, 8 and 9 and Tables 1 and 2 reveal shortcomings associated with the SFS 
representation of a flat ocean bottom. While both the ACS and SFS techniques reph- 
cated precisely the original contour, the latter introduced small errors in both the first 
and second-order derivatives. 

Figures 10, 11 and 12 and Tables 3 and 4 provide a comparison between the ACS 
and SFS representations of an ocean bottom which deepens (down slope) with constant 
slope. Table 3 and Figure 10 show that the ACS technique accurately modeled the ori- 
ginal function and its derivatives. However, Figure 10 and Table 4 highhght difficulties 
in the SFS representation of this same contour. This method accurately reproduced the 
original discrete bathymetric data values but introduced oscilhatory errors in the in- 
terpolated region between them. These small oscillations around the true function values 
were propagated through the subsequent derivative calculations generating large per- 
centage errors. 

Figures 13, 14 and 15 and Tables 5 and 6 provide an analogous comparison for an 
ocean bottom which shoals (up-slope) at a constant rate. The results observed for this 
model were identical to those for the constant down slope ocean bottom. 

Figures 16, 17 and 18 and Tables 7 and 8 pertain to the ACS and SFS represent- 
ations of the one period cosine ocean bottom. Both techniques accurately reproduced 
the fundamental function shape as shown by Figure 10 but the ACS approach had 
trouble with both the first and second-order derivatives. Figures 17 and 18 and Table 7 
refer. Additionally. both the ACS and SFS reconstructions of the first and second-order 
derivatives were characterized by large errors at the contour end points. 

Figures 19, 20 and 2] and Tables 9 and 10 provide a comparison between the ACS 
and SFS representations of the modified Witch of Agnesi ocean ( see [Ref. 6] ) bottom 
model. The SFS technique was incapable of reproducing this contour or its derivatives 
as shown in the referenced Figures and Table 10. Furthermore, the oscilliatory errors 
first observed with the constant up slope and down slope ocean bottoms were again 
present. The ACS approach yielded a far more accurate reproduction of the exact bot- 
tom function. However, both the first and second-order derivative values, although bet- 
ter than those reproduced bv SFS methods, were also unsatisfactory. 

Figures 22, 23 and 24 and Tables 11 and 12 relate to the ACS and SES represent- 
ations of the half catenary ( see [Ref. 6] ) ocean bottom. The problems observed with the 
SFS reconstruction of this contour were identical to those observed for the modified 
Witch of Agnesi bottom as seen in Table 12. The ACS method again provided a true 


reproduction of the exact mathematical bottom model. In addition, the first-order 
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Figure 7. 





8 9 8 8 8 R 9 


(W) 81 


Flat Ocean Bottom : Contour Reconstruction. 
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Flat Ocean Bottom : First-Order Derivatives. 
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: Second-Order Derivatives. 


Flat Ocean Bottom 


Figure 9. 
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Contour Reconstruction. 


Down-Slope Ocean Bottom : 


Figure 10. 
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Down-Slope Ocean Bottoni : First-Order Derivatives. 


Figure 11. 
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Up-Slope Ocean Dottom : Contour Reconstruction. 


Figure 13. 
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-Slope Ocean Dottom : First-Order Derivatives. 
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Figure 14. 


Ci 
єч 


530435 AHANI WELY 


1Ν{Ἴ45 “{8Π3 HUHH 














SNIHIOONS SOZONH' LNGHL IM 
Otala ONY 320439 SIN[Qd ULGO IJILIHINAS JQ 33840N [[ :SINIOd HI[IuO 0 4 ΗΝ UNIS IQ 
WOLLOS NJ3O0 Jd) Sd! LNMISNO) 4O JAC LYA(HOO GNZ 


(H1. 29NUX 


$a а 






`0 





(WIGA 40 NOLIU'QOdSJINI T1823ISM3WÜN YO4 SSNOINHISL 


Up-Slope Ocean Dottom : Second-Order Derivatives. 


Figure 15, 
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Ocean Bottom : Contour Reconstruction. 


Cosine 


Figure 16. 
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cean Bottom : First-Order Derivatives. 
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Figure 17, 
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Cosine Ocean Bottom : Second-Order Derivatives. 
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Figure 18. 


TECHNIQUES FOR NUMERICRL INTERPOLRTION OF YB(M) 


Figure 19. 
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Modified Witch Of Agnesi Ocean Bottom : Contour Reconstruction. 
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Modified Witch Of Agnesi Ocean Bottom : First-Order Derivatives. 


Figure 20. 
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Modified Witch Of Agnesi Ocean Bottom : Second-Order Derivatives. 
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CONSTANT 


AKIMA CUBIC SPLINE RESULTS 


DOWN-SLOPE BOTTOM 


Table 3. 
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(19VX3)1008A YOUN? Z 


00'025 
00'015 
00055 
00065 
00'09£ 
00'02$ 
ου ους 
00°O0S¢ 
00'0*£ 
00'0z$ 
00'02$ 
O00 ° OTE 
00°00¢ 
00062 
00'082 
00022 
00092 
00 052 
00°02 
00 #082 
00 022 


(H) 


( 113)8A 


00'025 
00°OT 
00°00 
00'06$ 
00'09$ 
00'02$ 
00'09$ 
00'0S$ 
00'0*$ 
00*'0££ 
00'02£ 
00 ‘078 
00°00¢ 
00'062 
00*' 0892 
00*'0422 
00'092 
00* 052 
00'0*52 
00 ‘082 
00 ‘022 


(H) 


(LIVXIIGA 


uo 
Qo o 


uou 
ЄС юк оо 
Q O QO O O Q QO QO QO QO O QO QO O O QO O O O r4 


ON O O O A O OA O O O 
NAMM Sst NNS 


00'0 


u) Ou) 
О r4 ri 


(ΜᾺ) 
JONVY 


114 3NITdS 9210402 VWINV : ЭПОІМНЭЗІ 


CONSTANT 


SERIES ERESULIS : 
DOWN-SLOPE BOTTOM 


FOURIER 


Table 4. 


HEH HIE 
X Xx XE XK € E ¢ 
E E E E € ¥ 3€ 3€ 
X X E E ¢ € € ¥ 
JE E ¢ IIE 
HEHEHE IEEE 
E E ЭЕ ЗЕ ЭЕ ЗЕ € 
XX XE EE ¥ ¥ 
X X3 X 3X 
E E E E E E E 
0000 

HEHE IE HHI 
E E HE IEE EIEN 
ЗЕ ЭЕ ЗЕ ЗЕ € ¢ 3€ 
JE FEI HIE HEH 
XX K X 3 
HH HHI HHH 
3¢ € € IE HEE 
X E € ¢ 3€ K€ 
E E € € 3€ E3 H 
X3X30000€ 


ЧОЧУ £ 


10-369 *0- 
10-2292 *0- 
20- 3094 0 
10-3497 '0 
20-382 *0- 
ТО- 3921 °0- 
20-3821'0 
T0-3901'0 
£0-3123'0- 
20-3646 '0- 
LT I900 O: 
20-3526 *0 
40-3725 0 
IL0-3901*0- 
20-3821'0- 
I0-3S21'0 
20-3£92'0 
το-349τ:0- 
20-3032 °0- 
T0-4292°0 
Τ0-36ας0 


(H/T) 


( 113)210Q08A 


00*3000'0 
00*3000'0 
00*3000'0 
00130000 
0030000 
00+ 3000 °0 
00*3000'0 
00+3000°0 
00+3000°0 
00*3000'0 
00*3000'0 
00* 3000'0 
00*3000'0 
00*3000'0 
00*3000'0 
00*3000'0 
00*3000'0 
00+3000°0 
00*3000'0 
00*3000'0 
00*3000'0 


( 12VX3)21008A 


II 


£91'994- 
992411 
gqI6'95t£ 
822 °S¢- 
CASE 
40° ST 
024 ° 8T 
597° 9- 
0S8'l9I- 
9461 
059 °SST 
9/6'l 
0s9g'l9I- 
59L°9- 
029° <8T 
5S50°ST 
24l'2t2- 
δαν 
ατό ους 
σαοσ᾽ατῖ 
£91'994- 


uou Z 


'SINTOd V1VG IWNISTYO dO ‘ON 


I0*3z2z1'0- 
00131250 
00538650 
00436210 
00+ 3592 '0- 
00*3022'0 
00134930 
00139910 
00135210- 
00139020 
00*3113'0 
00+3602°0 
00*3*»2I'0- 
00*t3981'0 
0013469530 
00130520 
00*3*592'0- 
00*3621'0 
00*398698'0 
00*3125'0 
L0*3£zl'0- 


( 114) 1008A 


00130920 
00130020 
00130020 
09130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 
00130020 


WOLLOG JJOISNMOG ANVISNOS 


000°0 
БЫА 
00070 
[GS o= 
000'0 
S65't 
000°0 
960'£- 
000°0 
0<0° < 
000°0 
922 `&- 
0000 
624't 
000'0 
692 b= 
000'0 
$50'4 
000'0 


g2 l- 


000°0 


(129VX3)100gA 80883 £ 


00°02 
29 00b 
00'00b 
6° 2 š 
00'03£ 
56285 
00.00 
91655 
00'052 
00'0*£ 
000565 
00°00¢ 
00'00£2 
5800ς 
00092 
20* 2S2 
00°092 
i? 292 
00 052 
25 °66ї 
00 ‘022 


(H) 


( 113)8A 


00:02 
00* 015 
00'005 
00'06£ 
00095 
00'02£ 
00'09£ 
00°0S¢ 
00'05£ 
00'0z£ 
00'02£ 
00'0T&£ 
00° 002 
00°062 
00°082 
00' 042 
00°092 
00°0S2 
00092 
00'0£2 
00'022 


(H) 


(LIVXIIAA 


00'I 
S6'0 
06'0 
αΌ 0 
090 
sZ°0 
040 
490 
090 
55 `0 
050 
3590 
050 
ας 0 
0£'0 
αὖ 0 
02°0 
Sl'0 
0l'0 
S0'0 
00°0 


CHM) 
39) Ἡ 


114 531135 W3lINDOJ 1VIl1VdS :3ñOINH331 
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CONSTANT 


CUBIC SPLINE RESULTS 


AKINIA 


Table 5. 


u0333 2 


UP-SLOPE BOTTOM. 


8T-3559°0 
8T-3822°0- 
LILI υ- 
81-39160 
11-3082 0 
££-32z21'0 
81-328 '0 
EI 0 
Ө1Т-3©©89°0 
££-32£1'0 
91-3£2£9'0 
t£- Jez Ὁ 
871-3258 '0 
&- 32571 ° 0 
81-328 '0 
££-32£l1'0 
981-3££29'0 
t£t£-Joel 0 
81-38 '0 
££-32£1'0 
87-38 '0 


(WT) 


( 114)210Q08A 


00+3000°0 
00*3000'0 
00*3000'0 
00*3000'0 
00*3000'0 
00*3000'0 
00+3000°0 
00+3000°0 
00+3000°0 
00+3000°0 
00+3000°0 
00+3000°0 
0030000 
00+ 3000 ° 0 
00+3000 °0 
00+ 3000 °0 
00*3000'0 
00*3000'0 
00*3000'0 
00*3000'0 
00*3000'0 


(W/T) 


(12VX3)210Q0dA  uONd3 Z 


TI 


00*3002'0- 
00*3002'0- 
00+3002 *0- 
00+3002°0- 
00+3002°0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00*3002"'0- 
00*3002"'0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00+3002°0- 
00+3002°0- 
00+3002°0- 
00*3002'0- 
00*3002'0- 
00+ 3002 *0- 
00+3002°0- 
00+3002°0- 


( 114)1O0QqA 


00*3002'0- 
00*3002'0- 
00*3002'0- 
00+3002°0- 
00*3002'0- 
00*3002'0- 
00* 3002'0- 
00*3002'0- 
00+3002°0- 
00+ 3002 °0- 
00+ 3002 °0- 
00+3002°0- 
00*3002'0- 
00*3002'0- 
00+3002°0- 
00*3002'0- 
00*3002'0- 
00+4002°0- 
00+3002°0- 
00+ 3002 *0- 
00+3002°0- 


( 199X3)10089A 


*SINIOd Vi¥d TVNI9I38O 40 “ΟΝ 


WO1108 3dOlSdü LINVLISNOO 


8403383 Z 


00*022 
00`0©2 
000597 
000532 
00°092 
000627 
00092 
00'062 
00°00¢ 
00°OTs 
00°02¢ 
00 0ες 
00°05© 
00°0S¢ 
00°09¢ 
00°0ZE 
00'09£ 
00'06€£ 
00'009 
00' 015 
00°02 


(H) 


να 


00' 022 
00'0£2 
00'0*2 
00' 082 
00' 092 
00'022 
00'082 
00'062 
00°00¢ 
00'0T€£ 
00065 
00 0τε 
00'05z 
00°0S¢ 
00°09¢ 
00°O0ZE 
00'08£ 
00°06¢ 
00°005 
00' 015 
00'02* 


(H) 


( 12VX3)98A 


(ИА) 
JINVY 


114 3ΝΙ145 этап VHIWV :3005INHO31 





πο ο p ишш αυ» 
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FOURIER SERIES RESULTS : CONSTANT UP-SLOPE 
BOTTOM. 


Table 6. 


эеэєэєэєэє эе Зее 
XXX X ¢ ¢ ¥ 3€ 
XX € ¢ ¢ X €3€ 
3X9 X0 
3¢ ¢ 3¢ 3¢ ¢ ¢ ¢ € 
HHH € € ¢ € 3€ 
є зе зе 3¢ X X X 
HH HH HH HH 
Hx AK 5€ ¢ € 
X3 FEE Xx HE HIE 
000'0 

HK KK RHE 
HK KKH KKH 
3¢ € € € € 5¢ ¢ € 
HA KI HHH 
339 969€ 
зе зе зе эе ¢ € ¢ € 
3¢ € 3¢ ¢ ¢ ¢ €3€ 
3¢ 3€ € € HEHE HH 
FEI ¢ € ¢ x3 
X300 X3 


ἩΟ8ὰ3 £ 


10-36ας”0 
T0-3292°0 
20-3092 *0- 
TO 9T 0> 
20-3282 °0 
το-3αστ0 
20-3821 '0- 
10-1901 0- 
©0-3125°0 
20-3626°0 
91-3922'0- 
20-3676 '0- 
£0=4025 0= 
I0-3901'0 
20-3821'0 
І0- 38921 °0- 
20-382 '0- 
Т0- 3291 ‘0 
20-30924 ° 0 
1023292 0= 
T0363: 0> 


(W/T) 


( LI3)21008A 


00+3000°0 
00*3000'0 
00+3000°0 
0030000 
00130000 
00130000 
0014300070 
004 30000 
00130000 
00130000 
00130000 
0030000 
00130000 
00+ 30000 
00+ 3000 *0 
00530000 
00130000 
00+ 3000 *0 
00330000 
00+ 3000 ° 0 
00+ 3000 °0 


( L2VX3)21008A 


tt 


ΣΉ] οσο 
S92 ат 
SI6'8^5z 
8L2°SE- 
CEE zam 
6S0°ST 
οσα °’ 8T 
ki 5- 
098° T9T- 
9461 
06S9°SST 
926°T 
05ο Τ9ἵ- 
992° 95 
02S8'£891 
23048 
24l'2£2- 
8/2 SF- 
616 ‘85 
S92°STT 
ESL 991- 


ЧОЧУЗ 7 


:*S1NIOd ViVGQ ПУМІЈІНО JO ‘ON 


I0*3££l'0 

00*3I£5'0- 
00* 39698'0- 
00+ 3621 °0- 
0041359920 

00*30£2'0- 
00+ 3299°0- 
00* 3981'0- 
00*3521'0 

00+ 3502 °0- 
00+3[IS°0- 
00*3502'0- 
00*35*521'0 

00+ 3981 °0- 
00+3/99*0- 
00+ 302 °0- 
00139920 

00+ 3621 '0- 
00+ 3868 '0- 
00*3T£5'0- 
TO+3ESET°O 


— ap a P a س‎ oe s s ae 


( LI43)1008A 


00+3002°0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00+3002°0- 
00*3002'0- 
00*3002'0- 
60*3002'0- 
00*3002 '0- 
00+3002 ° 0- 
00+3002 ' 0- 
00* 3002 ' 0- 
00*3002'0- 
00*3002'0- 
00*3002'0- 
00+3002°0- 
00+3002°0- 
00*3002'0- 


( L2VX3) LOQgA 


WOLLOG 3d01SdN LNVLONOD 


000°0 


SLC EL 


0000 
550'Z 
000'0 
692'*5- 
000°0 
ez E 
000'0 
9де = 
000°0 
0z0' £ 
000: 0 
9860'£- 
000'0 
56% '& 
0000 
οτι ο 
000 '0 
ἐνν 2 


κ 


00' 022 
£b6°66T 
000056 
Lo 292 
00°092 
10292 
00 '082 
98005 
00'00£ 
00°00 
00'02£ 
00°0%© 
00°05£ 
9T 672 
00'09£ 
£6'20£ 
00'08£z 
6£'24t£ 
00°00 
τα ΟΡ 
00'025 


(H) 


( LI4)8A 


00'022 
00'022 
00°02 
00°0S2 
00°092 
00'022 
00 ‘682 
00067 
00'00£ 
00'0I£ 
00'02£ 
00*0zz 
00°0%© 
00°O0S¢ 
00095 
00°0/© 
00'08£ 
00'06£ 
00°00$ 
00'0I5 
00' 025 


(W) 


( 1L2VX3)8A 


00'I 
S6'0 
060 
“80 


36 0 
04:0 
99 0 
09°0 
43° 0 
050 
q5'0 
05'0 
ας 0 
0z'0 
32'0 
02 °0 
ат`0 
От ‘о 
S0°0 
000 


(ΜᾺ) 
39NWVM 


414 511835 34314004 lVIIVdS :300INH231 
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ONE PERIOD 


τ: οᾱ- 
«040 
262 ° £9 
99£'92- 
220`°0% 
000'0 
13509 
7ης ις- 
£2l1'29- 
052'£1 
280792- 
052 ΣΤ 
950 ' 94 
952 °1t- 
13509 
0000 
14509 
99£'9c- 
BIS 95- 
«060 
22” 99- 


20-32% '0 
20-3992 '0 
Т0-32ТТ'0 
20-32 ° 0 
20-3258 °0 
21-3092 °0- 
20-3I16£'0- 
20-341£'0- 
20-32*2'0- 
20-3*58'0- 
20-35985'0- 
20-3658°0- 
I0-32II'0- 
20-3218 '0- 
20-316 '0- 
00430000 
20-316£'0 
20-32^*5£2'0 
20-32*£'0 
20-3952°0 
20-325 '0 


20-3062 *0 
20-3182 *0 
20-3629 '0 
20-35950 
20-362 ‘0 
Z1-389T 0- 
20-352 °0- 
20-3595°0- 
20-36£9'0- 
20-3192 *0- 
20-3062 °0- 
20-3154 ' 0- 
20-3629°0- 
20-3595 '0- 
20-3552'0- 
Z1-3522'0 
£0-3*5*52'0 
20-3595'0 
20-3629'0 
20-3192 *0 
20-3062 °0 


336 083833 Т0- 3622 °0- 


925 ОТ” 


632°2T 
eo [= 
5»£0'c2- 
656°0 
520°2- 
09g 2= 
É:9 SI 
¿bb° ¿- 
000°0 
¿bb ° ¿= 
££9'SI 
082'2- 
%0°2- 
6360 
5b£0'2- 
ДӨТ 
632°21 


beSCUILI- 
90090000 


00*325£'0- 
00*3628'0- 
1τ0:300Τ΄0- 
TO+3ZTT'0- 
TO#4ZeT*O- 
I0*3/ITl'0- 
00*3566'0- 
00*3*59'0- 
00+ 369 '0- 
51-3152 0 
00+ 369+ 0 
00+ 3598 '0 
00139660 
τ01341τ᾽0 
ї0+3/2Т`°0 
I0*3ZIl'0 
I0*3001'0 
00*3629'0 
00*3Z*5£'0 
I0-3624'0 


31-3999'0- 
00*399£2'0- 
00*36£/'0- 
I0*3201'0- 
I10*3021'0- 
I0*t3921'0- 
r0*302t'0- 
I0*320Tl'0- 
00*3624'0- 
00t399£'0- 
g31-35ST1'0 
00*399£'0 
0041366240 
TO+3201:0 
I0*30?T1'0 
I0*t3921'0 
I0*3021'0 
I0*3201'0 
00436640 
00+ 388+ °0 
00430000 


00°0S2 
995692 
02'992 
26 EE 
O02 99< 
00'055 
09g'IIS 
ΠΗ 
09'Il9 
83199 
00059 
85199 
089119 
PL 995 
οοτῖς 
ου 05υ 
02 839 
E6 E 
02'992 
99 2 
00°0S2 


00°0S2 
62'692 
02'992 
5b5'2£t 
06885 
0005» 
09'TIS 
99 293 
08 °ТТ9 
Тс '0%9 
00' 059 
I2'059 
091τ9 
ΠΗ 
09115 
00:05» 
06985 
DD’ CEE 
02'882 
67'692 
00'052 


00'I 
56 '0 
060 
«390 
090 
SZ°0 
02°90 
99 0 
09°0 
330 
030 
550 
0590 
St'0 
0£'0 
αό 0 
02 °` 0 
Sl'0 
ото 
S0'0 
00'0 
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AKIMA CUBIC SPLINE RESULTS 


COSINE OCEAN BOTTOM 


Table 7. 


(W/T) (H/T) (H) (H) (H9) 
(113)21008A ((12VX3)210089A 038893 Z (113)1009A (12VX3)1008A YONA Z (1143)9A (12VX3)0A JƏNVY 


τι :'SINIOd ViVQ TVNISTYO 40 ΟΝ 


114 3ΝΙ]ά5 314Π2 ΥΗΙΑΥ :“3ΠΒΙΝΗΩΘ3Ι. 


HOLIO8g 3AVM 3NISOO3 GOTY3d ЗНО 





ONE PERIOD COSINE 


FOURIER SERIES RESULTS 


OCEAN BOTTONI 


Table &. 


T99 95- 
SOL 2T 


σος ° 0T 
919 TIT 
062'9- 


HH Xx x3 xx 


Του ε- 
Tes z- 
9580 
I95'0 
SLS 0= 
T95°O 
35»9'0 
Ie3'2- 
[99 2- 


3€ 3€ €3€ t 


062'8- 
919'TI 
29£'01 


S07 2 [= 
Του ' 95- 


20-3262°0 
°0-3819°0 
20-3502°0 
ἐ0-38130 
20-3522'0 
<0- 3792 *0- 
20-3982 0- 
20-3295 '0- 
20-3559 °0- 
20-3537 '0- 
20-4982 °0- 
20-3592 *0- 
20-3509 '0- 
20-329%'0- 
20-3982 *0- 
£0-3132'0- 
20-3*22'0 
20-38135'0 
20-3502 ‘0 
20-3819 0 
20- 320+ ‘0 


( 114)21008dA 


20-3064'0 
20-3192 ‘0 
20-3690 
20-35900 
20-3552 '0 
21-3901 `°0- 
20-3552 ‘0- 
20-3595 '0- 
20-3629'0- 
20-3152 *0- 
20-3052 *0- 
20-3192 *0- 
20-369 °0- 
20-3590 °0- 
20-3552 *0- 
21-3622 '0 
20-35520 
20-3595 '0 
20-3659 '0 
20-37192 0 
20-3062 *0 


( 12VX3)21008A 


TI 


κκεκκκκ OOF ICICI 0- 


097 '9- 
022 79= 
ἴνα 0 
£90'I 
o91 g= 
985'0- 
301'0 
$2£'0 
56 I° 0- 
000°0 
%671°0- 
ασς 0 
“010 
885'0- 
o91 o> 
«401 
1540 
Ürc' v- 
0*51'9- 


HH IEE HH HIE 


uou 7 


:SAINIOd VIVO 1VNIS9INO JO ‘ON 


00*1^»92'0- 
00+3202°0- 
T0*320101'0- 
To+3tet’o- 
T0+392T*O- 
[0*36Tl'0- 
T0*3201'0- 
00*3152'0- 
00*39082z'0- 
9T=32¢¢ 0- 
00*398£'0 
00131540 
T0*3201'0 
I0*3611'0 
10534510 
Ι093Τ2Τ᾽0 
I0*3201'0 
00+3202°0 
00139950 
00*t3221'0 


( 114)10Q8A 


913998 0= 
00*398z2'0- 
00*3622'0- 
lo*3201'0- 
T0*302101'0- 
I0*3921'0- 
10+1021°0- 
IO+320[:0- 
00+36<2°0- 
00+3682°0- 
ST-3bST'0 
00+ 388 °0 
00536540 
ї0+4201°0 
1098305610 
I0*3921'0 
T0*t3021'0 
[T0*3201'0 
00136540 
00+ 388 ‘0 
00*3000'0 


(l12VX3)1008A  NOMNH3 Z 


WOLLO@ JAYM 3NISO2 GOIN3d ЗМО 


00 ‘082 
89 Το 
02 ‘882 
£8 TSE 
02°`60© 
62035 
ος τίς 
т әс 
08ττ9 
S2 059 
00'099 
S2° 059 
09 `II9 
z9 295 
03 [ITS 
6203» 
02835 
GEES 
06 986 
841927 
00°0S2 


(H) 


( LI418A 


00'092 
62°6S92 
06886 
bhb’ TT 
02'89£ 
00'08* 
ος τις 
997299 
08 °ТГ9 
Т2 ‘059 
00059 
T2°069 
08 *ТТ9 
947299 
οο τίς 
00'085 
0σ 885 
bose 
02882 
62'6832 
00°0S2 


( L2VX3)8A 


00'I 
36'0 
060 
39'0 
08'0 
32'0 
040 
α9 “0 
09:0 
33'0 
03'0 
9% '0 
090 
ας 0 
0z'0 
S2° 0 
020 
SI'0 
οτ’ο 
400 
00'0 


(WA) 
JONY Y 


114 531835 YJIXNOJ 1VILVdS :309INHO3L 
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WITCH OF 


AKIMA CUBIC SPLINE RESULTS 


AGNESI OCEAN BOTTOM 


Table 9. 


` 


CSL 
02£'4 
960'S2- 
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derivative results were accurate with a uniformly small percentage error as reflected in 
Table 11. Second-order derivative values generated by the ACS approach were unsatis- 
factory, although this is not immediately evident from Figure 24 because of the plot 
scale. 

Figures 25, 26 and 27 and Tables 13 and 14 give the ACS and SFS results for the 
parabolic ocean bottom model. The ACS method matched both the original function 
and its first and second-order derivatives precisely as shown in Table 13. The SFS rep- 
resentation was reasonably accurate although some oscillations were again in evidence 
at the contour end points. These oscillatrons account for the large boundary errors in the 
subsequent first and second-order derivative computations as reflected in Table 14 and 


illustrated бу Figures 26 and 27. 


B. NUMERICAL ENHANCEMENT OF THE SPATIAL FOURIER SERIES 
TECHNIQUE 

The results of the previous section clearly showed that the elementary SI'S approach 
was unsatisfactory. Modifications to the fundamental SIS technique were therefore in- 
vestigated and applied where appropriate to improve overall performance. 

The following numerical enhancements were studied for each of the original seven 
exact mathematical models. For purposes of brevitv, intermm results for the modified 
Witch of Agnesi contour only have been included, as this function reflected the greatest 
improvement, 

l. Lanczos Smoothing. 

The modified expression for the SFS coefficients given by Equation (2.11) and 
incorporating Lanczos smoothing was introduced into the SFS computer algorithm and 
the modeling process repeated. 

Figure 28 and Table 15 show that Lanczos smoothing was effective m the sup- 
pression of the between data point oscillations witnessed previously in Table 10 and 
Figure 19. The corresponding reduction in the magnitude of error for the first and sec- 
ond-order derivatives can be scen in Figures 29 and 30. One undesirable effect of this 
smoothing technique was the introduction of artificial curvatures at the end points of the 
bottom contour. The cause of this phenomenon is attributed to the apparent disconti- 
nuity of the function generated by the immediate commencement and cessation of data 
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Half Catenary Ocean Bottoim : Contour Reconstruction. 
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Parabolic Ocean Bottom : Contour Reconstruction. 


Figure 25. 
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Parabolic Ocean Bottum : First-Order Derivatives. 


Figure 26. 
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Derivatives: Fourier Series reconstruction with Lanczos smoothing. 
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Figure 30. Modified Witch Of Agnesi Ocean Bottom Second-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing. 
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2. Lanczos Smoothing With Synthetic Data. 

The ability of Lanczos smoothing to suppress oscillations generated bv the SFS 
technique gave it great utility. However, the artificial curvature at the end points of the 
contour was an undesirable effect which needed to be eliminated. It was unlikely that the 
problems associated with Lanczos smoothing could be countered by further modifying 
the expression for the SFS coefficients. Extension of the original bathymetric array by 
adding equally spaced synthetic data points was, however, viewed as a possible solution. 
The ratronale of this procedure was to force the undesirable features generated by 
Lanczos smoothing beyond the actual region of interest. 

Let each data point y;, be represented as a function of incremented range such 
that 


Ув = ув), {= 0,1,2,...(Х2ВРТ — 1) (930) 
with 
2,=iAz, *=0,1,2,...,(NZBPTS — 1) (3.2) 


where Az is the distance between each data point and NZBPTS ts the total number of 


Original discrete data. Now let 


Ув(21) = Ув(20) 


Hs == Az (52) 
апа 
Ja NZBPTS~1)) — Va wznprs-ay E 
a To πι Ад == аш: (3.4) 
Synthetic data points were then created by hnear extrapolation, that is, 
Ув(2_/) = ур(20) F n. و إ2‎ x 1,2 TS ТИР (079) 
where 
and 


¥al2vzaprs—149) 7 Jswzsprs-i) + Mr Zwzeprs—i4y +6 = 1,2,..., SINTHP. (3.7) 


95 


SYNTHP is the number of synthetic points to be generated. Figure 31 provides à 


diagramatic representation of this methodology. 


Range of 
Interest 


Synthetic Data Points 


Original Bathymetric Data 





Figure 31. Generation Of Synthetic Data: Synthetic duta points are generated 
by linear extrapolation. 


Figures 32, 33 and 34 and Table 16 reveal that the addition of five synthetic 
points to either end of the original data array was successful in shifting the end effects 
induced by Lanczos smoothing outside the horizontal range of interest. This resulted in 
a better curve fit for the modified Witch of Agnesi bottom and smaller errors in the 
subsequent first and second-order derivatives. One undesirable feature of the linear ex- 
trapolation method for simulating data was its tendency to suppress the natural curva- 
ture of the funcuon being modeled. This phenomenon is observed in Figure 32. 

Generating synthetic data by linear extrapolation was a valid technique when 
the nature of the function being modeled was relatively featureless. Помеуег, И маѕ пої 
suited to funcuons with substanual variauon in shape. Consequently, an alternative 
method involving the mirroring of the original data about the end points of the contour 


Was created. This inethod is described by 


Yal — 2) =yp(2) im 12,-,SYNTIHIP (3.8) 
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TRANGE IKM) 
NUMBER OF SYNMHETIC DATA POINTS BCTORC AND AFTER: S 


TECHNIQUES FOR NUMERICAL INTERPOLATION OF YBI(M) 


MODIFIED WITCH OF AGNES) OCEAN BOTTOM 


ORIGINAL NUMBER OF DATA POINTS: 2) 


WITH LANTZ0S SMOOTHING 





Figure 32. Modified Witch Of Agnesi Ocean — Bottom — Contour 
Reconstruction: Pourier Series reconstruction with Lanczos smooth- 


ing and lincar extrapolation. 
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Fieure 33. Modified Witch Of Agnesi Ocean Bottom First-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 
and linear extrapolation. 
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Figure 34. Modified Witch Of Agnesi Ocean Bottom Second-Order 
Derivatives: [Fourier Series reconstruction with Lanczos smoothing 
and linear extrapolation. 
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and 
Ув (хавртѕ-1+0)) = Ув мевртѕ-1-)) » 1i l2, 5 EEINTHTE. (3.9) 


Figure 35 illustrates this mirroring process. 


Range of 
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Synthetic Data Points 


O 
al 


@ = Original Bathymetric Data 


Figure 35. Generation of Synthetic Data: Synthetic data is generated by mirror 
maging the original data points. 


l'igures 36, 37 and 38 and Table 17 give the results for the SES modeling of the 
modified Witch Of Agnes: contour using synthetic data points derived from the miror- 
Ing technique. Further improvement in the overall curve fit and accuracy of the first and 
second-order derivatives was observed. However, some leveling of the bottom contour 
curvature was sull prevalent. The cause of this phenomenon was attributed to the ap- 
parent function disconununty arising from the mirroring process. This problem 15 identi- 
fied in Figure 35 where nurror naging at the right contour boundary has created a spike 
in the data. In an effort to compensate for this, the SIS technique has smoothed this 
feature, thereby suppressing the true funcuon curvature at the end points of the range 
oi Milcrest. 

To resolve this problem of curvature suppression, the nurror imaging technique 


given by Equations (3.8) and (3.9) was altered as follows: 
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Figure 36. Modified Witch Of Agnesi Ocean Bottom Contour 
Reconstruction: Fourier Series reconstruction with Lanczos smooth- 
ing and mirrored data. 
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Figure 37. Modified Witch Of Agnesi Ocean Bottom First-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 
and nurrored data. 
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Figure 38. Modified Witch of Agnesi Ocean Bottom Second-Order 
Derivatives: Fourier Series reconstruction with Lanczos sraoothing 
and mirrored data. 
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Ув( = 2) = ув(да-р), {= 1,2,...5ҮМГНР (3.10) 
and 
YelZvzeprs-i+)) =YalZwzeprs—p) ἐπ 1,2,...SYNTHP. (3.11) 
Figure 39 illustrates this modification. 
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Figure 39. Generation of Synthetic Data: Synthetic data is generated from a 
shifted mirror imaging of the original data. 


Figures 40, 41 and 42 and Table 18 reflect only a marginal unprovement in the 
modified Witch of Agnesi curve fit as a consequence of using the inodified mirror imag- 
ing procedure. Additionally, for other exact mathematical models such as the one period 
cosine contour, this procedure proved detrimental, again because of curvature sup- 
pression, 

An alternative means of eliminating function discontinuities created by mirror 
imaging was investigated. This problem was most pronounced when the slope between 


successive data points increased toward the end points of the contour. A new method 
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Figure 40. Modified Witch of Agnesi Ocean Bottom Contour 
Reconstruction: Fourier Series reconstruction with Lanczos smooth- 
ing and shifted mirrored data. 
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Figure 4l. Modified Witch Of Agnesi Ocean Bottom First-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 
and shifted nurrored data. 
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SPATIAL FOURIER SERIES 
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Figure 42. Modified Witch of Agnesi Ocean Bottom Second-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 
and shifted murrored data. 
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SHIFTED MIRROR 
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IMAGED DATA. 
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for creating synthetic data was therefore developed to take advantage of this change in 


contour slope. Let 


= yg(21) – Ув(со) 


Mio = ДЕ (3.12) 
121 = 2800001) , (3.13) 
Ув(2(үғвртУ-1)) — ya(2wzeprs-2)) 
Dini eme рә Aa (3.14) 
and 
T т (3.15) 
where again Az 1s the distance between each data point. If 
mo st (3.16) 
and, or 
[пл < | rt, | ; (3.17) 
then synthetic data 1s generated at the appropriate end of the profile as follows: 
ув(2_1) = yp(2o) * GAz(mig — iA) , i= 122,.,8Y NTHP (3.18) 


and, or 
Js UqwzaPrs-142) 7 JE Zwznprs-1) * GAz(mg, — iàm,)) , i— 1,2,.,8 YNTHP (3.19) 
where 
Ат = туу = ht) (3.20) 
апа 


Amp = Mg) — Me. (3.21) 
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A comparison of the end points in Figures 35 and 43 demonstrate the utility of this new 


method in preserving curvature and eliminating apparent function discontinuities. 
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synthetic Data Points 


Original Bathymetric Data 





Figure 43. Generation of Synthetic Data: Synthetic data is generated by revers- 
ing the rate of change of slope of the original boundary data. 


Enhancing the original bathymetric data with synthetic points proved beneficial 
to the SFS modeling of the bottom contour, but no one procedure for creating this ar- 
tificial data dominated. Each of the methods investigated showed merit under specific 
circumstances. Consequently, the SES computer algorithm was modified to incorporate 


the three methods outlined previously. If 
[Ан = [от 1 < (3.22) 
and, or 
lAmgl 9 Img, — mal «e, (3.23) 
where € is some arbitrarily small number, then the synthetic data was constructed by 


linear extrapolation as given in Equations (3.5) and (3.7). However, if 


72 


zao | mns, (3.24) 
and, or 
Е < | ту, | 5 (5,25) 


then mirror imaging as given by Equations (3.8) and (3.9) is used to generate the data. 
Finally, if the conditions given by Equations (3.16) or (3.17) hold then the method em- 
ploving the reversal of the rate of change of slope is activated to simulate data at the 


appropriate end point. 


C. SPATIAL FOURIER SERIES ANALYSIS OF THE AKIMA CUBIC SPLINE 
REPRESENTATION OF THE OCEAN BOTTOM CONTOUR. 

Results from the SFS technique contained significant errors in the reconstruction 
of both the original contour and first and second-order derivatives despite Lanczos 
smoothing and the optimization of methods for generating synthetic data. However, for 
the more complex ocean bottom models such as the modified Witch of Agnesi and one 
period cosine functions, SFS methods yielded the more accurate derivative values. 

The ACS method provided precise reproductions of the original test functions, but 
the subsequent first and second-order derivative values Were reliable only for the flat, 
constant slope and parabolic ocean bottoms. Efforts were therefore made to incorporate 
the desired aspects of both methods. 

The ACS representation of the ocean bottom: model was discretized at ten times the 
sampling frequency used to generate the original bathyinetric data. Synthetic data was 
then constructed as described in Section B of Chapter 3, with the exception that Az was 
now based upon this new sampling frequency. For exainple, where previously 5 synthetic 
data points covered a 500 m range, now 50 synthetic points describe this same range. A 
SFS representation of the entire contour was then obtained. This approach was applied 
to each of the seven exact mathematical functions. 

Figures 44, 45 and 46 and Table 19 provide results for the SFS representation of the 
flat ocean bottom generated by the ACS method. The original contour is replicated 
precisely and the small errors observed in both the first and second-order derivatives 
approach the limits of machine precision. However, while these errors are very small, 
they are actually larger than those generated by the fundamental SFS reconstruction of 


the flat bottom without Lanczos smoothing or synthetic data. Table 2 refers. 
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Figure 44. Flat Ocean Bottom Contour Reconstruction: Spatial Fourier Series 
representation of the Akima Cubic Spline flat bottom model. 
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Figure 45. Flat Ocean Bottom First-Order Derivatives: Spatial Fourier Series 
representauon of the Akima Cubic Spline flat bottom model. 
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Figure 46. Flat Ocean Bottom Second-Order Derivatives: Spatual Fourier Series 
representation of the Akima Cubic Spline flat bottom model. 
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Figures 47, 48 and 49 and Tuble 20 give the results for the SFS representation of the 
constant down-slope ocean bottom generated by the ACS technique. The percentage 
errors recorded for both the bottom contour and the first and second-order derivatives 
are two orders of magnitude lower than the initial results given at Table 4. This re- 
duction, although significant, is not sufficient to eliminate residual errors likely to ad- 
versely impact on ray acoustic computations. 

Figures 50, 51 and 52 and Table 21 provide results for the SIS representation of the 
constant up-slope ocean bottom model. The observations made for this case are identical 
to those recorded for the constant down-slope bottom. 

Figures 53, 54 and 55 and Table 22 give results for the SFS representation of the 
one period cosine ocean bottom. The errors for both the contour reconstruction and its 
first and second-order derivatives were larger than those recorded in the earhest test 
provided in Table 8. Furthermore, the sampling of the ACS model for this bottom shape 
not only failed to assist SFS in more accurately reproducing the contour, but caused it 
to replicate anomolies in the ACS reconstruction of the first and second-order deriva- 
tives. 

Figures 56, 57 and 58 and Table 23 provide results for the SI'S representation of the 
modified Witch of Agnesi ocean bottom. This contour was used to demonstrate the de- 
velopment of successive numerical enhancements to the fundamental SI'S technique, and 
therefore, gives the best opportunity for a complete comparative analysis. Significant 
improvement Was observed in the modeling of this function with SFS methods employ- 
ing Lanczos smoothing and synthetic data created by nurror imaging. A comparison of 
the results in Tables 17 and 23 show that further reductions in error were achieved by 
the SFS approach to modeling the ACS representation of the modified Witch of Agnesi 
bottoni. However the recorded inprovements in the first and second-order derivatives 
remain insuflicient for the purposes of existing ray acoustics algorithms. 

Figures 59, 60 and 61 and Table 24 give results for the SI'S representation of the 
half catenary ocean bottom. Significant reductions in the percentage errors of the con- 
tour were achieved. For the first and second-order derivatives reductions in percentage 
error of three orders of magnitude were common. However, the residual errors remained 
sizeable and prohibited the use of these results in further computations. 

Figures 62, 63 and 64 and Table 25 provide results for the SI'S representation of the 
parabolic ocean bottom. The method of using the ACS model of the test function as 
input to the SFS procedure again proved beneficial with significant reductions in the 


magnitude of the percentage errors being observed for the contour. Reproductions of the 
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Figure 47. Down-Slope Ocean Bottom Contour Reconstruction: Spatial Fourier 


Series reconstruction of the Akima Cubic Spline down-slope bottom 
model. 
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Figure 48. Down-Slope Ocean Dottoin First-Order Derivatives: Spatial Fourier 
Series representation of the Akima Cubic Spline down-slope bottom 
model. 
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Figure 49. Down-Slope Ocean Bottom Second-Order Derivatives: Spaual Fourier 
Series representation of the Akima Cubic Spline down-slope bottom 
model. 
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Figure 50. Up-Slope Ocean Bottom Contour Reconstruction: Spatial Fourier Se- 
ries representation of the Akima Cubic Spline up-slope bottom model. 
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Figure 531. Up-Slope Ocean Bottom First-Order Derivatives: Spaual l'ourier Se- 
rics representation of the Akima Cubic Spline up-slope bottom model. 
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Figure 32. Up-Sleupe Ocean Bottom Second-Order Derivatives: Spatial Fourier 
Series representation of the Akima Cubic Spline up-slope bottom 
model. 
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Figure 53. One Period Cosine Ocean Bottom Contour Reconstruction: Spatial 
Fourier Series representation of the Akima Cubic Spline cosine bottom 
model. 
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Figure 54. One Period Cosine Ocean Bottom First-Order Derivatives: Spatial 
Fourier Series representation of the Akima Cubic Spline cosine bottom 
model. 
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Figure 55. One Period Cosine Ocean Bottom Second-Order Derivatives: Spatial 
Fouricr Series representation of the Akima Cubic Spline cosine bottom 
model. 
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Figure 56. Modified Witch Of Agnesi Ocean Bottom Contour 
Reconstruction: Spatial Fourier Series representation of the Akima 
Cubic Spline Witch of Agnesi bottom model. 


88 


SERIES ` 


-- 


- 
— —À o 


SPATIAL FOURTER 





5 
A 
С 
ο 
8 
E 
Ё 


=оз = = = = = = = 


0.6 0.8 


0.7 


SPATIAL FOURICR SORICS F171 OBTAINED BY SAMPLING PRIMA CUBIC SPLINE FIT WITH SAMPLING RATE 10 


NUMBER OF SYNTHETIC DATA POINTS BEFORE AND НЕТЕР: 50 


TECHNIQUES FOR NUMERICAL INTERPOLATION OF YIR) 





Ë 
О ОО з e Ἢ 
E 
IE 
° É 
= 
Ы 
NT E e 
Vm 
B 
= 
| το 
JE 
b B 
M umm... p c m 
£ > Ñ 
> ` ° 
ЕЕ 
"σα 
, 985 
8 i £ р h š b = £ 
|a aea T mc | 
в ë ë ë ë б с οὗ ë 6 e 
Un 


Figure 57. Modified Witch of Agnesi Ocean Bottom First-Order 
Derivatives: Spatial Fourier Series representation of the Akima Cubic 
Spline Witch of Agnesi bottom model. 
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Figure 58. Modified Witch of Agnesi Ocean Bottom Second-Order 
Derivatives: Spatial Fourier Series representation of the Akima Cubic 
Spline Witch of Agnesi bottom model. 
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Figure 59. Half Catenary Ocean Bottom Contour Reconstruction: Spatial Fourier 
Series representation of the Akimna Cubic Spline half catenary bottom. 
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first and second-order derivatives Were very accurate over the major portion of the range 
of interest, but were characterized by large errors at the outer end points of this region, 
as illustrated in Figure 64. Tlus phenomenon is attributed to the reversal of slope niethod 


for generating svnthetic data. 


D. ALTERNATIVE INTERIM MODELING TECHNIQUES 

Efforts to model arbitrarily shaped ocean bottom contours by SFS methods yielded 
unsatisfactory results in every case studied despite three numerical modifications de- 
signed to enhance the performance of this technique. ACS results for the contours were 
accurate, but the subsequent derivative calculations were erroneous. The need for a 
program which could accurately model both the bottom and its derivatives remained 
paramount. Without this, the precise calculation of reflection angles for acoustic rays 
impinging an undulating sea floor could not be realized. 

Motivated by the above requirements, attempts were made to salvage the ACS 
technique by utilizing its accurate contour reproduction and applying alternate numer- 
ical methods to obtain reliable values for the first and second-order derivatives. The 
following central differencing algorithms, extracted froin [Ref. 7], were utilized: 


| — ув(лу + 2А) -- 8rg(zg + А) — Syg(zg — h) + Vp (29 — 2h) 
ro В\20 В\20 T p Zo B\ZO (3.26) 


and 


— Ypo + 2A) + 16vglzo + h) — 30rg(zg) -- 16rp(zg — А) — ypg(zg — 2A) 


E (3.27) 
Zh 


У в(20) = 
Where A is some arbitrarily small number, nominally 1/100 of the average spacing 
between each original data point. 

The accuracy of Equations (3.26) and (3.27) was verified in separate tests using se- 
veral mathematical functions. These relationships were then applied to the ACS recon- 
structions of the seven exact mathematical ocean bottom models listed in Section 1 of 
this chapter. The resultant first and second-order derivatives were an improvement on 
previous methods, but below the precision expected. These poor results are attributed 
to the nature of the ACS program and the manner in which the snterpolating splines 
between each data point are generated. 

Having abandoned both the SES and ACS imethods for modeling the ocean bottom, 


attention was focused upon a computer routine which provided a divided differences 
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Figure 62. Parabolic Ocean Bottom Contour Reconstruction: Spatial Fourier Se- 
ries representation of the Akima Cubic Spline parabolic bottom. 


100 


SPATIAL FOURICR SERIES FIT OBTRINCD BY SAMPLING RKIMA CUBIC SPLINE FIT WITH SAMPLING RATE 10 


- 
> 
L 2 
E E 
£ Е 
= z 
E : 
= 
© 
= z 
LJ © 
Ξ 
E E 
: : 
E š = 
= š 
L.) e 
= 5 
Е а, 
= š 
tL Б 
Ë 
О = 
= ° š 
5 T 
E ΠΕ 
5.5 ع‎ 
Ἵν 
E25 
ba > 
c d x 
VES 





Figure 63. Parabolic Ocean Bottom First-Order Derivatives: Spatial Fourier Se- 
rics representation of the Akima Cubic Spline parabolic bottom. 
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interpolating polynonual fit to unordered, randonily spaced data. This program , written 
bv Gerald and Wheatley [Ref. 7], was tested on the seven exact mathematical bottom 
models. Accurate curve fits were achieved; however, more care was required in sampling 
each contour to avoid artificial curvatures in the fit which are often characteristic of this 
numerical technique. The subsequent first and second-order derivatives generated by 
applying Equations (3.26) and (3.27) to the polynomial fit were precise. Therefore, these 
equations and the polvnomial coefficients program formed the basis of an interim mod- 
eling technique which would enable verification of the reflection angle algorithm for an 


arbitrarily shaped ocean bottom. 
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IV. ACOUSTIC RAY REFLECTIONS FROM ARBITRARILY SHAPED 
BOTTOM CONTOURS 


The reflection angle algorithms given by Equations (2.13), (2.14), (2.15) and (2.35) 
were incorporated into an existing acoustic ray trace program. This Fortran program 
had previously been limited to the flat bottom case. Seven exact mathematical functions 
where then used to model ocean bathymetry. One acoustic ray was launched for each 
bottom contour and the computer generated angles of reflection and bottom slope 
compared with the corresponding theoretical values. An isospeed ocean was assumed in 
each case In order to interpret the results easily. 

The seven exact mathematical functions employed were: 

e [lat Bottom (у,(2) = 1000 т) 

@ Constant down-slope bottom (y,(z) = 600.0 + 0.04z m) 

e Constant up-slope bottom (y,(z) = 1000.0 — 0.04z m) 

e One period cosine bottom (y,(z) = 250.0 — 50.0 cos(0.0002z7z) m) 


Bodied Witch of Agnesi bottom (y,(2) = 200.0 —-——1 =) _ m) 
: 0.000122 + 2500.0 


e Half catenary bottom (y,(z) = 315.0 — 30.0(¢7/6000.0 + @-2/6000.0) туу) 


(z — 5000.0)? 
75000.0 


Figures 65 through 7$ illustrate the ocean bottom models and the propagated 


e Parabolic bottom (у,(2) = 250.0 + m). 


acoustic ray for each of the test cases. A comparison between theory and the computer 
generated results for each model are given in Tables 26 through 32 inclusively. For each 
table, the terms y^ and £^ represent the exact values for the angle of bottom slope and 
acoustic ray reflection angle, respectively. The terms y? and е represent the estimates 
generated by the ray trace program and reflection angle algorithm for these same pa- 
rameters. The bottom and surface bounces are denoted by BB and SB, respectively. 
These results were uniformly accurate with almost all of the recorded errors for the 
computed reflection angles, [2 below 0.1%. Errors in the estimated bottom slope, ў°, 
were also reasonably small thus verifying the central differencing methods for determin- 
ing the first-order derivatives. The parabolic ocean bottom model yielded the worst re- 
sults with errors in the reflection angle, €°, averaging 0.08%. The modified Witch of 


Agnesi and half catenary models generated the largest errors in angle of slope, y?, but 
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the magnitudes of these angles were sufficiently small such that little impact on the 
computed angles of reflection resulted. 

Arithmetic confirmation of the reflection angle algorithin is not readily determined 
for a ray launched in an ocean medium with speed of sound varying as a function of 
depth and for an arbitrarily shaped ocean bottom. However, two such cases were inves- 
tigated as a means of visual verification of the effects of a ray incident upon an arbi- 
trarily shaped ocean bottom. The one period cosine and modified Witch of Agnesi 
models were utilized after shght modification to linut the range of depth. Figures 79 and 


80 refer. The sound-speed profile is given by Table 33. 
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Figure 65. — Flat Ocean Bottom: This profile was used to test the reflection angle 
algorithm for an isospeed ocean. 
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Figure 66. Acoustic Ray Impinging a Flat Bottom: This figure shows an acoustic 
ray launched at depth 15 m with initial angle £, = 45° and propagating 
in an isospeed ocean with a flat bottom. 
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Figure 67. Constant Down-Slope Bottom: This profile was used to test the re- 
flection angle algorithin for an isospeed ocean. 
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Figure 68. Acoustic Ray impinging a Constant Down-Slope Bottom: ‘This figure 
shows an acoustic ray launched at depth 15 m with initial angle 
Во = 20* and propagating in an isospeed ocean with a constant down- 
slope bottom. 


Table 26. REFLECTION ANGLES: FLAT BOTTOM WITH /, = 45° 























Table 27. REFLECTION ANGLES : DOWN-SLOPE BOTTOM WITH f, — 20* 


rop ШШЕ 


[ 6BB asm | 2291 | 229 | 000 | 45.97 [45197 | 9.00 — 
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Figure 69. Constant Up-Slope Bottom: This profile was used to test the reflection 
angle algorithm for an isospeed ocean. 
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Figure 70. Acoustic Ray Impinging a Constant Up-Slope Bottom: This figure 
shows an acoustic ray launched at depth 15 m with initial angle 
fi, = 55° and propagating in an 1sospeed ocean with a constant up-slope 
bottom. 
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Table 28. REFLECTION ANGLES : UP-SLOPE BOTIOM WITH 6, = 55° 


BB | 3315 | 2291 | 2291 | 0.00 [| 48428 | 38428 | 000 — 
SB [aa | 0000 | 0.00 | 0.00) 45.838 | 2353535 | 0.00 
ET 
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Figure 71. Опе Period Cosine Ocean Bottom: This profile was used to test the 
reflection angle algonthm for an 1sospeed ocean. 
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Figure 72 Acoustic Ray Impinging a Cosine Bottom: This figure shows an 
acoustic ray launched at depth 15 m with inital angle Û, = 30° and 
propagating in an isospeed ocean. 


Table 29. REFLECTION ANGLES : COSINE BOTTOM WITH £, = 30° 


ISB 0224 [0.000 | 0.000 | 0.00 | 30.242 | 30.248 | -0.019 
7288 | 0.262 | 0.000 | 0.000 | 0.00 [31.007 | 31.009 | 0906 — 
[ssp ons | 0.000 [0.000 | 0.00 | 32.298 | 32.297 | 0903 — 
[app | osa | 0910 | 0910 | 0.00 | 33.208 | 33.207 | -0.003 
"asp | ess | oo | 0.000 | 0.00 | 34418 | 34418 [ 0900 — 
L— sp | 2023 | oo | 0.000 | 0.00 | 42366 | 41367 | 0002 _ 
вв | 2766 | 1774 | 1776 | on [arsos | 47893 | 6000 
osn | 3070 | 0000 | 0000 | 000 | 49:67 | 49:68 | 0002 — 
απ | 5399 | 0.000 | 0000 | 0.00 $2605 | 54605 | 0.000 — 
Lis ss озю | oo | 0.000 | 0.00 | 52849 | 52850 | 0002 - 
[ase | 6596 | 0000 | 0.000 | 0.00 | 49351 | 493554 | 0906 — 
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Figure 73. Modified Witch Of Agnesi Ocean Bottom: This profile was used to 
test the reflection angle algorithm for an isospeed ocean. 
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Figure 74. Acoustic Ray Impinging a Modified Witch Of Agnesi Bottom: This 
figure shows an acoustic ray launched at depth 15 m with inititial angle 
6, = 70° and propagating in an isospeed ocean with a Witch of Agnes 
bottom. 
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Table 30. REFLECTION ANGLES : WITCH OF AGNESI BOTTON WITH 
Po = 70° 


πμ [το | 0.000 | 0.000 | 000 — 70.644 | 70.663 | 0.027 — 
0.336 | -030 | 72439 
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Pictire-75. “Halt Catenary-Otern_ Bottom: Eis or one Wis seu toate st tie 
flection angle alyorithm for an isospeed occan. 
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Figure 76. Acoustic Ray Impinging a Half Catenary Bottom: This figure shows 
an acoustic ray launched at depth 15 m with initial angle Û, = 60° and 
propagating in an isospeed ocean with a half catenary bottom. 


Table 31. REFLECTION ANGLES : HALF CATENARY BOTTOM WITH 
fa = 60° 


ise ова | 0000 | 0000_| 000 | 39321 | 59.931 | 0017 
Ces | 1294 | 0425 | 0125 | 0o00 | 39.800 | 59.805 | 0.008 
OBB | 2457 | 0210 | 0214 | 150 — [59461 | 59466 | 0908 — 
Ciee | 7622 | osa | osa | 0.00 [30.287 | 50256 | -0.002 
Cse [ass | 0.000 [0.000 [000] 49.46] 49345 | 0002 - 
Cisse [8.675 | oo | 0.000 | 0o00 [45072 | 45.073 | 0902 — 
Lass | 9.015 | 0.000 | 0.000 | 0.00 [42698 200 | 0905 — 
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Figure 77. Parabolic Ocean Bottom: This profile was used to test the reflection 
angle algorithm for an 1sospeed ocean. 
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Figure 78. Acoustic Ray Impinging a Parabolic Bottom: ‘This figure shows an 
acoustic ray launched at depth 15 m with initial angle £, = 70° and 
propagating In an isospeed ocean with a parabohe bottom. 


Table 32. REFLECTION ANGLES : PARABOLIC BOTTOM WITH /, = 70° 


Range Ὁ Λο 0 ° ^o T 


298 | sas | 0.000 | 090 | 909 — [61255 | st313 | 9078 — 
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Figure 79. Acoustic Ray Impinging A Cosine Bottom: This figure shows an 
acoustic ray launched at depth 5 m with initial angle f, = 85° and 
propagating In an ocean medium where sound speed varies as à func- 
tion of depth. 
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Figure 80. Acoustic Ray linpinging A Modified Witch ОГ Agnesi Bottum: ‘This 
figure shows an acoustic ray launched at depth 5 m with initial angle 
Въ = 85° апа propagating in an ocean medium where speed of sound is 
a function of depth. 
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Table 33. SOUND SPEED PROFILE 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The IMSL Version 10 Akima Cubic Spline (ACS) computer computer program had 
proven useful in modeling ocean bottom contours, but was unable to provide accurate 
first and second-order derivatives for the one period cosine and modified Witch of 
Agnesi ocean bottom models. The application of Spatial Fourier Series (SES) methods 
as an alternative means of modeling the one-dimensional ocean floor was considered a 
viable solution to this problem. However, despite several numerical enhancements to the 
fundamental technique, including Lanczos smoothing and artificial data simulation to 
eliminate end point anomalies, Spatial Fourier Series failed to achieve the accuracy re- 
quired for subsequent acoustic ray and reflection angle computations. 

Despite the aforementioned difficulties, Spatial Fourier Series methods may vet 
prove to be a viable numerical tool in modeling ocean bathymetries, but its immediate 
utility, based upon the analvsis conducted, 1s doubtful. 

The divided differences polynomial curve fitting Fortran routine yielded results 
which were equal to, or better than, both the ACS and SFS modeling techniques that it 
replaced. Similarly, the central differencing method for determining the first and 
second-order derivatives of the ocean bottom provided the most accurate estimates when 
used in conjunction with this alternate approach to bottom contour modeling. 
Notwithstanding the foregoing observations, this polynomial routine also had deficien- 
cies, the most notable being the introduction of artificial oscillations generated following 
an imprudent choice of discrete bathymetric data to describe the contour. 

The reflection angle algorithm developed in Chapter 2 for an arbitrarily shaped 
one-dimensional, ocean bottom were validated by the computer simulations given in 
Chapters 

One of the realizations of this research was that no one numerical modeling method 
studied was definitive. Further investigation of bottom contour modeling techniques, 
such as the method of least squares, 1s warranted. This work represents precursory ef- 
forts in the development of a two-dimensional ocean bottom model which would enable 
more advanced study of ray acoustic theory in an ocean medium where sound-spced 


varies as a function of all three spatial coordinates. 
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